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Abstract. 


An electromagnetic pulse (EHP) model is developed which allows a 
quick computation of the time development of the electric fields qene- 


rated by a high altitude nuclear burst. The model is based on the 


MN eee men ic, Sin 


Karzas-Latter high frequency approximation for high altitude EMP, which 


Cescribes fields generated by Compton electrons interacting with the 


The proper choice of a gamma time output funee 


earth's magnetic field. 


tion, which can be Integrated in closed form, and a smal} angle approxi- 


mation, made in the expressions for the Compton currents and alr conduc- 


tivity, eliminate the time consuming numerical Integrations usually 


necessary in EMP models to compute the Compton currents and air con~ 


ductivity. This results in a considerable savings tn computation time. 


The mode’ is presented in a manner which is simple to use but sttl! 


allows the vari tton of the major theoretical parameters in the prob- 


lem. 


A. simplified model of electron collision frequency as a function 


of electric field strength is given which enables the model to predict 


accurate results for nuclear weapon gamma ylelds up to at least 100 Kt. 


The results predicted by this EMP model compare to within 5.5% with 


results from the Air Force Weapons Laboratory CHEMP computer code. 


The computation time using the presented model on a CDC 6600 compu- 


ter Is typically 5 sec or less for a 5 shake computation period in 


steps of .1 shake. 
The model presented should be useful for both classroom Instruc- 


tion and nuclear vulnerability/survivability studies and analysis prob< 
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(1. Theory 


Cverview 


Since the solution to the EMP problem is actually the solution of 
a classical electromagnetic theory problem, the derivation of the moJel 
equations reduces to putting Naxwell's equations into 3 convenient form. 
This is essentially accomplished by expressing Maxwell's equations in 
spherical coordinates and transforming to a retarded time frame. One 
q must also develop expressions for the currents and conductivities of 
4 


4 the system in the absorption region. The general equations d.reribing . 


2 
se ees. 


the high altitude model of Karzas and Latter have been derived in great 


; detail by Chapman (Ref 3). Only the major points of the derivation will : . 


be given here. The system origin is assumed to be at the burst point 
with detonation at time t = 0. This geometry is illustrated in Fig. 1. : a: 

The key points to be remembered in this model are: 

1. Each gamma ray gives rise to one downward traveling 

Compton electron. 4 

B 2. The electrons are turned by the earth's magnetic Fleld s a 
. giving rise to a centrifugal acceleration, q 3 
3. The relativistic electrons radiate eneray in their 


forward direction. 
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4, The gamma rays and the EMP radiation travel at the q 
a same speed. This leads to constructive interference 


of the radiation from each of the electrons. 
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a Particle Densities 


;: The gamma rays from cne nuclear weapon travel in a straight line 


to a point where they produce Compton electrons. At any given point r 
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the aumber of ganmas which interact to produce Compton electrons fs 


2 
exact 208s 
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a g(r) = - - 2 Q1) 
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where A(r) is the mean free path of gasma rays to produce Compton elec- 


trons, Y is the gamna yield of the weapon in electron volts (eV), and E 


is the mean gamma energy in eV, 


a. Equation (1) may also be called the radia) distribution function, 


or an attenuation function for interacting gamva rays. 


The : term is 


The 4ar2 


term accounts for the divergence of the gamma rays as the radius r ts 


the total number of gamma rays available from the weapon. 


increased while the remaining terms account for the reduction in gammas 
due to thir absorption in the atmosphere, based on the mean free path. 


It is assumed that the gama mean free path varies as the exponen- 


eats 


tial atmosphere. 


This gives the functional relationship between A and 
r: 


A(r) = Ag exp [(HOB - Fr cos A)sos (2) 
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where 


se 


tae 


49 = gamma mean free path at standard pressure 
HOB = height of burst in km above the earth's surface 


r = radial distance from the burst point to the point of interest 


PAN) BYERS 


> 
A = angle between the position vector r and the vertical 


S$ = atmospheric scale height 


With this assumption, Eq (1) can be integrated and tecomes 


re 1 (. Ss _ HOBN FF ricos AN _ 
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(3) 
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Now if f(t) is the time distribution function of the weapon yield, 


the rate of Compton electrons, ne, produced at a given ooint r and time 


t is given by 


mee alr) F(e-£) (4) 


Each Compton electron produces through inelastic scattering events 


several secondary electrons which form the basis for the conductivity 


of the atmo~’"ere. As in the Karzas-Latter approach, each Compton elec- 


tron is assumed to have a constant speed, Vo, throughout the ranae, R, 


of the electron which is a function of altitude. This allows the life- 


time to be expressed as R/NVo. If each Compton electron produces secon- 


dary electrons at a constant rate, the rate of secondary electron, ns, 


production is 


dng E¢/33 eV 4 
q a «RNS (5) 


where Ec is the energy of the Compton electron and 33 eV ts the average 
ionization energy per air molecule (Ref 2). 
Considering the differential current produced by the Compton elec- 


trons, it can be shown (Ref 3) that the Compton current and the number 


areal 


of Compton electrons are given by 


No 


+ R - 
JS(t) =-eglr) f — V(r*) f ‘ - 7+ AG >) dt’ (6) 
FS . 


c 


R/V, - 
nels) = ofr) f° (r= 7 # SED) a (7) 
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where 


- 


t the time since the creation of the Compton electron 
X(t’) = the radial distance the Compton electron has traveled 
e = the magnitude of the electron charge 
The quantity +t is aenerally known as retarded time. 
it then follows from Eqs (4) and (7) that the number of secondary 


electrons is 


Vv t RIV 
ng(t) = = g(r) f f S o¢ G- -1t7' + 
=o 0 


where q is E,/33 eV. 


Currents and Conductivity 

tn the Karzas-Latter theory, the speed of the Compton electrons ts 
considered to be a constant, however, there is an acceleration cue to 
the seomaanetic field. The general equation of motion for an electron 


tn this case is 


d > > + > > 
" Wa-e (E+V X B) - mv 


electron rest mass 
electron velocity 
electric field 
magnetic field 
the electron collision frequency 


(1 = (Vg/e)2)~% 
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If the relativistic motion of the electron is considered, only the 
VX B term is imnortant and in soberical coordinates, the expressions 


for the veloctty components become (Ref 3) 
Vp = Vo (sin? 9 cos ut + cos’ 0) (19) 
Vg 4 Vy (cos @ sin @ cos ut - sin 8 cos 6) (1) 
Vg = Vo (sin 9 sin ut) (12) 
where w i- the cyclotron frequency for an electron and is given by 


(13) 


with B, the maanitude of the geomagnetic field. 


From Eq (10), X(r°) is found to be 
- # sin wt” ri 
X(t ) = Vo (sin? o SRE + x cos? 0) 


The Compton currents may now be written as 


iY 


RIV 
JE (t) = - eg(r) Vo f ° [Ee (T) (cos? 6 + sin? 9 cos wt)]} dt* (15) 
° 


¢ i RNo 
J (t) = - eg(r) Vy J {fF (T) stn 6 cos 6 (cos wr* - 1)J dr” 
t) 


(16) 


¢ RN 
ae =~ eg(r) Vy {F(T) sin 9 sin or’) dr? (17) 


where 


x 
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In a Similar manner Eq (8) becomes 


a6 


t R/V, 
ng(z) = Me g(r) J [f° eer) are] are 
7 ° 


where 


v v ee 
Teer C1 - 2 cos? @) 177 + sin? 9 SHOT 
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wo 


(20) 


Equations (35), (16), (17), (18), and (20) may be simplified if 
the factor ur Is assumed to be small. In this case a Taylor series 
expansion for the sin and cos of wt, including only first and second 
order terms, is 


sin wt = wT (21) 


wer? 
cos wt = 1 + (22) 


expressions for the Compton currents now become 


R/V, 2 RIV 
- ea(r) Vo f s 2 f(T) dr? - sin? 6 mae | © Ze (T) dr’ 
Lo 


1°26 (T) dt’ 


R/Vo 
° 


2 
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where ~ = Vo | 
¢ 
In a like manner Eq (20) becomes 
To at’ > (1 - 8) 17” (27) 


An expression for the conductivity may alse be found by using the 
equation of sotion for the secondary electrons. These electrons are in 
the thermal reaions with energies ranaing from about 10-15 eV to the 
ampient energy. {t should be remembered for later use that the ambient 
energy of the secondary electrons is dependent on the electric field 
present. For consideration here, it is assumed that y % 1 and also that 
the change of velocity with time is smail conpared to the other terms 
in Ea (9) so that fi may be neglected. Also, with low velocities, the 
¥ x z term is small compared to the remaining terms and “ay also be 
neglected. Then the velocity of the secondary electrons is 


Va-fue 
MWe 


Using Eq (28), the current due to the secondary electrons ts 


>, > 2 > 
JS(t) = - e Ving(t) = ae E ng(t) (29) 


>. > 
Comparing Eq (29) to J5 = o£, an expression for the conductivity 


o(t) = £— ng(t) (30) 
bed of 


Equations (19), (23), (24), (25), (26). (27), and (30) provide the 


destred expressions for the Compton currents and the conductivity. 


pi ler ate a Sides = alam Noreen age Z pat or 6c , on 


_ LEER Retest carereccceengs cs 


ray 


L£ 


GEP/PH/75-33 


sae 
wher 


Se eee ln a wena lla cera Ul 


Field Equat tons 


Maxwell's equations is rationalized MKS units are 


VXEo - 38 (31) 


(32) 


(33) 


(34) 


where q, is the total charge density and J is the total current density. 


'n addition to these equations the continuity of charge requires that 


- aay ~> 
—+ 93 2 0 
a 7 (35) 


> > 
Combining these equations to separate £ and B and transforming 


them into spherical coordinates and into the retarded time frame (Ref 


+ > 
3) the relations for E and 8 become 


-~ 2 4 >> 7 > 
- VE + Up —— Ved + = MQ 
C£g fv 


(36) 


a. 72133 * oo - 
tor ora fre) + uy (J - up J) | 


tn the Karzas-Latter model, only the time derivative portion of 


Eqs (36) and (37) are kept since the current variation with distance 


is slow compared to the variation in time for the high frequency con- 
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ponents. Also the fields and currents vary rapidly in time. This 


approximation is valid for about 100 shakes. 


Th $s same high freauency 


approximation is used here. In addition, the radial component of the 


field is dropped since it is weak compared to the transverse components 


and contributes only a very tow frequency siqnal (Ref 2). The equa- 


tions for the transverse components are 


(rEog) * up Yea] 7° (38) 4 
3 213 Pes) 2 
- [ eae (re,) = Jy ] 0 (39) 


(40) 


The currents in Eqs (38), (39), and (40) are total currents. The 


total currents are given by 


c 
Jog * Jot o(r) Eg (41) 


Substitution of Eq (41) into Eqs (32), (39), and (40) and inteora- 


tion over time gives 
c 
+ Uodg + to a(t) Eg = 0 


+ Body + uo a(t) Eg = 0 


u " 
+ 258 + Lat) fy = 0 
c c 
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Equations (42) and (43) are tn a fora which can be solved. The 
terme needed for solution of these equations are the air conductivity 
and the transverse components of the Conpton currents. The air con- 
ductivity may be found by using Eqs (19) and (30). The transverse com- 


ponents of the Compton currents may be found by using £qs (24) and (25). 
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TYPICAL LONG PULSE 
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~ 2. Tyoical Long Gama Output Pulse 
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TYPICAL NARROW PULSE 
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Fiq. 3. Tyoical Narrow Gama Output Pulse 


SLRS 2a 
A A Ace nel Tele m o Tpe Ro an antl ene ne a ne ge 


GEP /PU/ 75-13 


WV, Results. 


To test the validity of the HAEMP model, the computed values of 


this model were compared to available equivalent model values from the 


AFWL code for EMP calculations knowr as CHEMP (Ref 4). 


these calculations was: 


The basic set of conditions used for 


u 


target location around zero 


height of burst = 100 Km a 


3 aauss or 3 (10)75 wh/m2 


rs 


geomagnetic field 


0.6 degree 


inclination angle 


Compton electron recoil eneray = .75 MeV 


The aamma yield was va ‘ed and the number of steps taken for the 


numerical integration was varied according to the gamna yield, with 


Other parameters varied for 


more steps taken for the higher yields. 


examination of peak field values were burst height, aecmagnetic fiela, 


A preionization level was also considered. 


and pulse shape. 
The available data fron the CHEMP(N) code, which is CHEMP run with 


non-sel f-consistent calculations, was computed using a pulse of the 


form of Eq (46) with a = 1079 and & = 107 and the sane geometry as 


given above. This pulse shane is almost identical to that of Eq (48) 


Usina this 


See Fig. 2 on page VE. 


with a = 107 and 8 = 3.7 (10)8. 


pulse shape, a range of gamaa yields from .01 Kt to 100 Kt was used to 


The peak field values are plotted in 


calculate the EMP field vatues. 


The values taken from the CHEMP(N) code are annotated by X 


Fig. 4. 


marks. These peak field values show a maximum difference of 5.5% 


around the | Kt case, and a difference of less than 2% for al! other 


known cases. 
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Fig. 6. Peak Electric Field as a Function of Burst Height 
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HEIGHT GF BURST = 100 KM 
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Fig. 7. Peak Electric Fields for Tio Gamma Pulse Shanes 
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HEIGHT GF BURST = 100 KM 
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Fig. 8. Effect of a .03 Kt Ganma Yield Precursor Furst 
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Fig. 9. Effect of Geomaanetic Field Streagth on Peak Fields 
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